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The intrinsic dynamics of a system with open decay channels is described by an effective non- 
Hermitian Hamiltonian which at the same time allows one to find the external dynamics, - reaction 
cross sections. We discuss ways of incorporating this approach into the shell model context. Several 
examples of increasing complexity, from schematic models to realistic nuclear calculations (chain of 
oxygen isotopes), are presented. The approach is capable of describing a multitude of phenomena 
in a unified way combining physics of structure and reactions. Self-consistency of calculations and 
threshold energy dependence of the coupling to the continuum are crucial for the description of 
loosely bound states. 
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I. INTRODUCTION 


The center of interest in modern nuclear physics has 
recently moved toward nuclei far from the valley of stabil¬ 
ity. Weakly bound nuclei cannot be fully described in the 
limited framework of the shell model with a discrete en¬ 
ergy spectrum. Even the properties of their bound states 
reflect the proximity of the continuum. Loosely bound 
nucleons create an extended spatial structure that de¬ 
termines the results of possible reactions so that nearly 
all excitation mechanisms break up the nucleus. The 
standard approaches of many-body theory, such as the 
Hartree-Fock-Bogoliubov mean field and random phase 
approximation, necessarily include virtual and real exci¬ 
tations to the continuum. The Borromean cases of ®He, 
®Be, and ^^Li, when the system can be considered to be 
made of three clusters with all two-body subsystems be¬ 
ing unbound, are very sensitive to the continuum physics. 
This is the area where the conventional division of nu¬ 
clear physics into “structure” and “reactions” becomes 
inappropriate, and the two views of the process, from 
the inside (structure and properties of bound states) and 
from the outside (cross sections of reactions), should be 
recombined. 

The broad success of the nuclear shell model with ef¬ 
fective interactions urges one to look for ways to incorpo¬ 
rate the rich experience accumulated in the shell model 
into a more general context which would properly include 
the continuum part. We will not discuss below the most 
complicated task in this direction, namely the problem 
of the effective interaction. It is virtually unknown what 
should be an effective interaction of quasiparticles in the 
restricted shell model space which includes the contin¬ 
uum. The work that started with the Brueckner-Bethe 
theory of the G-matrix has to be reviewed and modified. 

Our goal here is much less ambitious. We would like to 
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demonstrate the new qualitative effects that emerge with 
the simple reformulation of the shell model in terms of 
an effective non-Hermitian and energy-dependent Hamil¬ 
tonian describing the “inside” view of the dynamics in a 
many-body system of interacting particles coupled to and 
through the decay channels. For our limited purpose be¬ 
low we assume that the effective interaction of the shell 
model can be simply readjusted to the new problem, al¬ 
though in fact it can be non-Hermitian by itself. 


The description with the aid of an effective non- 
Hermitian Hamiltonian is well known going back to the 
classical Weisskopf-Wigner damping theory |]^, works in 
atomic physics by Rice 0 and Fano 0 and projection 
formalism by Feshbach |^. The consistent formulation 
of the approach was given in the book by Mahaux and 
Weidenmiiller in application to processes with one 
particle in the continuum. This gave rise to the shell 
model embedded in the continuum iJ recently revived 
[§ ^ 1^ for the description of loosely bound nuclei. An¬ 
other direction of development was related to the descrip¬ 
tion of statistical and chaotic phenomena in nuclear reac- 
tion s |ll[ and generalization of random matrix theory 
[]I^ |14|| for unstable systems. The detailed study of the 
effective non-Hermitian Hamiltonian revealed new col¬ 
lective phenomena II, [T^ with bright manifestations in 
nuclear physics of low |lf JT^jand intermediate ||^ 
energies, atomic physics] 2C)|, |21|, molecular physics 
quantum chemistry |^, and condensed matter physics 
(p^ , 1^ . The basic origin of this collectivity is the 

same as in the Dicke superradiance [^, coherent cou¬ 
pling of intrinsic states through common decay channels 
(common radiation field of atoms confined to a small vol¬ 
ume in the Dicke case). The ideas related to this ap¬ 
proach were used for an analysis of experimental data, 
especially in two-state examples taken from nu¬ 

clear and mesonic physics as well as from the microwave 
cavity experiments ^ ^ . Here we follow a generic 
path of the shell model, adding continuum effects by in- 
cluding explicitly non-Hermitian terms in the Hamilto¬ 
nian M, H. 
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II. NON-HERMITIAN HAMILTONIAN 

We will not repeat here the full derivation of the effec¬ 
tive non-Hermitian Hamiltonian that can be achieved by 
separating the full Hilbert space into the intrinsic part 
and the continua and eliminating the continuum part 
with the aid of projection operators. This procedure was 
addressed in detail by many authors, see for example 
10, |3^. We label intrinsic states by 1,2,..., and the con¬ 
tinuum channels by a, b, c,.... The matrix elements of the 
effective intrinsic Hamiltonian can be written as 

'^12 = ^^12 + Ai 2 — - Wi 2 , ( 1 ) 

where H is an internal, let us say, a standard shell-model 
part, and the last two terms, which in general are func¬ 
tions of running total energy E, are generated by the 
exclusion of the continuum. 

The imaginary part W {E) originates from the real pro¬ 
cesses of decay to channels that are open at a given en¬ 
ergy. It is represented by the residues of the on-shell 
terms corresponding to the delta-functions coming from 
the energy conservation and causality requirement im¬ 
posed on the energy denominators, E E^~^'>. The 
quantity W has a factorized form, 

^ HJAf, ( 2 ) 

c;open 


amplitudes A‘1 as those determining the entrance and exit 
channels in Eq. (^). This guarantees the unitarity of the 
^-matrix since the virtual processes of evolution of the 
open system to and from the continuum channels are 
included in all orders in the propagator. Indeed, if one 
introduces the intrinsic Hermitian propagator {E — EI)~^ 
and the second order Hermitian scattering amplitude. 
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where the propagation does not include the coupling to 
continuum, the full scattering amplitude is given by the 
geometric series (the hats mark the operators in the chan¬ 
nel space). 


T = 


K 


1 -h {il2)K 

Then the scattering matrix is explicitly unitary, 


S = 


^1/2 l__(z/2)A 

l + {i/2)k 


( 6 ) 


(7) 


The diagonalization of the non-Hermitian Hamiltonian 
d^) produces the complex eigenvalues 


Ea,{E) = E^{E) - -T^{E), 


( 8 ) 


where the decay amplitudes A\{E) are the matrix ele¬ 
ments of the original total Hermitian Hamiltonian be¬ 
tween the states |I) and |c; E) of different subspaces; the 
normalization coefficients are included in the definition 
of A\. The second term of Eq. (|^, Ai 2 (£’), originates 
from the principal value of the same expression and cor¬ 
responds to the virtual off-shell processes taking place 
via the continuum. Therefore it includes contributions 
from all, open and closed, channels. For the system in¬ 
variant under time reversal, one can use a real intrinsic 
basis, where the matrix elements H 12 , A 12 and A\ can be 
taken real. 

The same effective Hamiltonian (]^) determines the 
scattering amplitude and the reaction cross sections. The 
relation between the inside and outside views was stud¬ 
ied in [|| 1^, 1^, ^ 1^, Q . The scattering matrix in the 
channel space describing the b a process is given by 

sab ^ (ga)l/ 2 (^ab _ T^b^(^^by/2 ^ ( 3 ) 


= ^—I (4) 

^ 1 \E-nJ,, ^ ^ ’ 


where the real, Ea, and imaginary, Tq, parts are func¬ 
tions of running real energy E. Without explicit energy 
dependence of the effective Hamiltonian, these eigenval¬ 
ues would provide the unstable states with a pure expo¬ 
nential decay law oc exp(—Tat). The presence of energy 
dependence violates the exponential decay, and the ac¬ 
tual quasistationary states are found at real energies Ea 
determined by the self-consistency condition 

Ea{Ea) = Ea. (9) 

The line-shape is not Breit-Wigner but we still call 
I'aiEa) the width of the resonance a. In what follows 
we omit the tilde sign for Ea if it does not lead to a 
confusion. 

In the region of interest, namely for loosely bound sys¬ 
tems, the main energy dependence comes from the prox¬ 
imity of thresholds as was stressed in Refs. 0,11 • The 
channel c is open only if the total energy E is above the 
threshold energy for this channel. The decay ampli¬ 
tudes associated with the channel c contain therefore the 
step factor Q{E — E^‘^'>) and can be written as 

Al=al{E)Q{E-E^‘^'>), (10) 


Here s“ = exp(2*(5o) stands for the smooth scattering 
phase coming from remote resonances not accounted for 
explicitly. The propagator {E — Ti)~^ in the scattering 
amplitude T“^ does not depend on a specific reaction and 
contains the full effective Hamiltonian (|^) with the same 


where af (E) is a smooth function of energy that falls off 
to zero when energy decreases to the threshold value. For 
a single-particle decay channel, it can be parameterized 
[ 0 , 01 as proportional to the square root of the pene¬ 
trability in this channel. 
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The real part A of the effective potential can be written 
as the principal value integral 

= - E / (11) 

71 - Jem E - E' 

Under the same assumption of a non-singular character 
of of, the matrix elements ( 0 ) also have a smooth energy 
dependence with no singularities near threshold, and can 
be approximated by energy-independent quantities as in 

0ii- 

One formal conclusion concerning the existence of 
bound and unbound states can be reached just from the 
way the theory is constructed. If the conventional shell 
model with a purely discrete spectrum (no coupling to 
the continuum) predicts a state with energy below all 
decay thresholds, this state will remain bound in the full 
calculation with the decay amplitudes included. Indeed, 
all widths depend on the total energy and vanish below 
thresholds so that the old solution is still valid. However, 
this statement is formal since it assumes that the reac¬ 
tion thresholds are known beforehand. In fact, they have 
to be determined consistently for the chain of nuclides 
relevant to the reactions under consideration. 

Certainly, there are limitations in the applicability of 
the effective non-Hermitian Hamiltonian method in the 
form outlined in the present article. As energy increases, 
a rapid growth of a number of interfering open channels 
makes this approach impractical. The approximations of 
a different type can be then introduced in the general 
framework, see for example that directly lead 

to quantum kinetics of statistical reactions . We also 
deliberately limit ourselves here by taking into account 
only the energy-dependence associated with threshold 
and resonance phenomena although the smooth “poten¬ 
tial” scattering part could be included without significant 
difficulties via the entrance and exit scattering phases 
hidden in the factors s“ and of Eq. (||). The full 
energy dependence was discussed, in particular, in Refs. 

1^ . The main physical assumption made here is that 
the states under consideration are close to threshold and, 
at relatively low energy, only few open channels are re¬ 
ally essential. The deviations resulting from violations of 
these conditions have been studied numerically in [ p4[ . 
For the purpose of this paper, namely for the develop¬ 
ment of shell model methods intended for the description 
of low-lying states in the nuclear systems near the border 
of stability, the presumed conditions are usually fulhlled 
to within a sufficient accuracy. 


III. SHELL MODEL APPROXIMATION 


We view the Hamiltonian (|^) as a sum of three terms. 


n = H° + V -^W , 


( 12 ) 


where we assume that the intrinsic Hermitian part EI°+V 
consists of independent particle energies, 

H° = '^e^ala^, (13) 

and the effective Hermitian interaction V. As a renor¬ 
malization of the standard shell model interaction, the 
Hermitian matrix elements A 12 , Eq. (0), generated by 
the virtual coupling through continuum, can be incorpo¬ 
rated into the operator V. The approximation of energy 
independence of the operator A and, as a result, energy 
independence of V , also used in previous works ^ , 
can be easily removed. 

In order to formulate physical problems in the spirit 
close to the conventional shell model, we start from the 
basis states |$), the eigenstates of H°, that are Slater de¬ 
terminants of the TO-scheme or their linear combinations 
projected onto correct values of total spin J and isospin 
T. For the next step we consider the “unperturbed” part 
of the Hamiltonian which includes the independent par¬ 
ticle part H° and the imaginary part —{i/2)W 

n° = H°- ^W. (14) 

The diagonalization of the Hamiltonian (|lj) along with 
the self-consistent solution of Eq. gives new eigen- 
vectors |$) either with complex energies (||) or as stable 
configurations on the real ener^ axis. The factorized 
nature of the operator W, Eq. (|2|), that is preserved by 
orthogonal transformations, and the presence of symme¬ 
tries may bring additional simplifications. In some spe¬ 
cial cases, 'H° remains diagonal in the original basis |<I>). 
In these situations the meaning of the amplitudes is the 
most clear, being related to the single-particle decay into 
continuum. 

Thus, for an isolated single-particle level \i^) embedded 
in the continuum, the unperturbed real energy is Ecore + 
e^. If the only open channel, c => is associated with the 
emission of the particle v, the imaginary part W leads to 
the width 


7. = KP (15) 

for any configuration which consists of the particle on 
the level Iv) and an arbitrary state of the stable core 
(no interaction between them at this stage). Threshold 
energy is determined by the core configuration. Simi¬ 
larly, in the case of several single-particle levels ly embed¬ 
ded in the continuum, the single-particle decay channels, 
opened for a specihc conhguration |<I>) with occupation 
numbers ni/(4>) = 0 or 1 , result in the width 

7 (^) = 

V 

If there are several single-particle levels jt'i) with the 
same exact quantum numbers j'^T, the situation is more 
complicated and in general E[° and cannot be simul¬ 
taneously diagonalized. Apart from the particle emission 
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from a given single-particle state we have now also the in¬ 
teraction through continuum given by the off-diagonal el¬ 
ements of Wi 2 ■ Here the Dicke collectivization and redis¬ 
tribution of the widths [@, |^, are possible, see below, 
and already at this stage, with no residual shell-model 
interaction V, we need to diagonalize the non-Hermitian 
Hamiltonian . 

This situation is almost certainly present in cases 
where two-particle emission is possible from different ini¬ 
tial configurations. For example, a zero spin pair can be 
emitted from a few j-levels leading to a final state with 
the same quantum numbers and therefore into the same 
decay channel. Here the coupling through the continuum 
may be very important. 

In the most general case, because of the energy depen¬ 
dence in the amplitudes Af, the positions of the resulting 
quasistationary eigenstates of 1-1° in the complex plane 
should be determined avoiding false solutions that can 
appear due to the possible non-analytic energy depen¬ 
dence at thresholds. One can take only those complex 
poles that can be traced back to the real axis (indepen¬ 
dent particle states 1$)) in the case of the closed chan¬ 
nels. Another new feature is that, generally speaking, the 
threshold energies are not known a-priori. They are to 
be calculated self-consistently comparing total energies 
of the parent and daughter nuclei taken in the same ap¬ 
proximation. But this must be done only on the next step 
when the residual shell-model interaction is accounted 
for. 

Finally, we include an effective interaction V that in¬ 
troduces the mixing of (in general unstable) shell model 
configurations 1$). As a result, instead of original inde¬ 
pendent particle states, we obtain the states I'k) which 
characterize both the intrinsic structure and possible de¬ 
cay channels in the fully interacting system. Since the 
Hamiltonian is non-Hermitian, the resonance “energies” 
(^ move in the complex plane relative to their initial po¬ 
sitions given by eigenvalues of H° (states |l>)) and this 
dynamics may be quite complicated driving some states 
back to stability. 

Another important “intermediate” Hamiltonian If° -|- 
V describes the case when all decay widths are “switched 
off” and obviously corresponds to the standard shell 
model. Below we denote shell-model many-body eigen¬ 
states as I'I's.ni.). Of course, in practice it is not necessary 
to make a two-step diagonalization, and the intermedi¬ 
ate steps with the wave functions |$) or |'ks.m.) can be 
avoided. 

It is known, see for instance [^, that the eigenstates 
of a non-Hermitian Hamiltonian form a biorthogonal sys¬ 
tem which would allow one to study the observable char¬ 
acteristics of unstable states along with the reaction cross 
sections and dynamics transformed to the time domain. 
Below we show several examples ordered by increase of 
complexity, from very schematic to more realistic. The 
selected cases illustrate the diversity of the physical phe¬ 
nomena that can appear in unstable many-body systems 
and can be described by the method of an effective non- 



bound orbitals 

FIG. 1: 

Hermitian Hamiltonian. As mentioned in the Introduc¬ 
tion, in this exploratory study we assume the effective in¬ 
teraction to be known and its matrix elements V 12 (taken 
in a basis of stable states and including A 12 ) to be real. 

IV. SINGLE-PARTICLE DECAY IN A 
MANY-BODY SYSTEM 

A. One single-particle level in the continuum; 
energy-independent continuum width 

We start with the simplest problem (a similar example 
was shown in Ref. |Q ). Consider a set of single-particle 
energies with the upper of them lying above the par¬ 
ticle emission threshold see Fig. |^, where a system 
of Q single-particle levels is presented with for 

ly = — 1 and for u = ^. For simplic¬ 

ity we assume here that the levels are equidistant on the 
real axis and not degenerate; later we add the Kramers 
double degeneracy. We assume that all single-particle 
emission channels have different quantum numbers and 
cannot be coupled through continuum. Thus, we have 
only one initial non-zero single-particle width expressed 
with the aid of the complex single-particle energy 

'I 

Now we use this set of single-particle levels as the basis 
for forming the many-body configurations |$) as Slater 
determinants with all possible distributions of the occu¬ 
pancies rii,. Finally we switch on a real two-body inter¬ 
action V. 

In order to characterize the generic results which are 
insensitive to specific peculiarities of the residual interac¬ 
tion, in this example we use a system of fV = 4 fermions 
with the equidistant spectrum e,y of H = 8 orbitals and 
random (Gaussian distributed) matrix elements of the 
two-body interaction. We solve this problem by diago¬ 
nalizing the complex matrix of the full many-body Hamil¬ 
tonian H. The results are shown in Fig. || as the dynam¬ 
ics of the complex eigenvalues evolving as a function of 
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the only variable parameter, the single-particle width 7 , 
taken here as an energy-independent number. 

At 7 = 0 (a normal shell-model limit), the many-body 
states |d>s.m.) obtained with the real residual interaction 
are stable, and their spectrum can be represented by 
the points on the real axis. As 7 increases, all states 
|d^s.m.) 1 ^') acquire widths and move into the com¬ 

plex plane. This means that, because of the unrestricted 
configuration mixing, any many-body eigenstate |'I') con¬ 
tains an admixture of configurations with the occupied 
upper orbital, and therefore it is capable of decay. In the 
limit of small 7 , in agreement with usual perturbation 
theory it is expected that Idt) « Idts.m.). Thus, the decay 
width of a many-body state T is, similarly to the case 
of Eq. dil), determined by the spectroscopic factor of a 
progenitor stable state I'I's.m.) that is given as an occu¬ 
pation probability of a decaying single-particle orbital 


r(4-) =7n,(«'s.n..). (18) 

This natural picture breaks down once the value of 7 
becomes comparable to the level spacing D along the 
real axis, and the internal dynamics gets affected by the 
continuum in a non-perturbative way. 

In this example 7 , and therefore the effective Hamilto¬ 
nian 7i, is independent of running energy E . This makes 
the trace of a conserved quantity resulting in 

The real part of the trace remains con¬ 

stant and is independent of 7 . In Eq. (H) we counted 
the number of non-interacting configurations available for 
[N —1) particles if the decaying orbital is occupied. Con¬ 
servation of the trace of Ti^ and the fact that the imagi¬ 
nary part W is diagonal make ( 7 )“^ E(d>)r('I') also 

a 7 -independent constant. 

In the 7 —> 0 limit, the occupancy of the decaying or¬ 
bital governs the distribution of widths in Eq. ( p^ . This 
can be generalized by introducing the dynamic occupa¬ 
tion numbers 




dm; 7 ) 

dj 


( 20 ) 


These parameters describe how at a given 7 an infinites¬ 
imal increase of the initial single-particle width 7 is dis¬ 
tributed among the many-body states dt. According to 
Eq. ( 0 ), no(^', 7 ) is independent of 7 . Despite all 

the resemblance to occupation numbers the numbers (|^ ) 
can be negative. One can also introduce generalized spec¬ 
troscopic factors r('I'; 7)/7 that are always positive and 
bound between 0 and 1 , but being cumulative quantities 
they would be less sensitive to dynamical features. 

With further increase of 7 the picture in Fig. looks 
paradoxical |Q. The many-body states are clearly di¬ 
vided into two groups. The complex energies of the first 


group rapidly move away from the real axis revealing 
large widths. At the same time the states of the sec¬ 
ond group turn back to the real axis and keep only tiny 
widths, whence becoming long-lived. Similar phenom¬ 
ena are known for a long time from various schematic 
studies 1^, ^ and versions of the shell model with the 
continuum effects included, for example |^, Q. 

The puzzle is readily resolved since with 7 increasing 
we come to the situation where the imaginary term dom¬ 
inates the dynamics and therefore classifies the eigen¬ 
states by their relation to decay rather than by real en¬ 
ergy. Any superposition of configurations with a consid¬ 
erable amplitude of the occupied unstable level 8 under¬ 
goes fast decay. Such states constitute the first group. 
The superpositions of the configurations with the empty 
level 8 become eigenstates of the second group and corre¬ 
spond to long-lived compound states. It is easy to calcu¬ 
late the dimensions of the two groups. The total number 
of states Af for 4 fermions on 8 non-degenerate orbitals is 
8!/(4!)^ = 70. The first group contains the states where 
the level 8 is occupied and the remaining 3 particles are 
distributed over 7 stable levels; the corresponding dimen¬ 
sion is 7!/(3!4!) = 35, so that the eigenstates are divided 
in this case evenly between the two groups, in agreement 
with Fig. ||. 

We see that the strong coupling of intrinsic states to 
the continuum produces a natural segregation of pro¬ 
cesses into fast direct reactions and slow compound nu¬ 
cleus reactions. Here the direct processes are of single¬ 
particle nature, and at large 7 the corresponding width 
of each state of the first group is close to 7 . Fig. ||, 
where the segregation of the effective occupation num¬ 
bers is shown as a function of 7 , confirms the division 
of states into short-lived with a fully occupied decaying 
single-particle orbital, ng = 1 , and compound ones where 
ng = 0. The analog of this segregation effect has been 
experimentally seen in microwave cavities, although 
for realistic systems the situation is more complicated, 
mainly because of the uncertainty in the effective num¬ 
ber of open channels ^ . 

In order to further quantify this phase transition we 
introduce a parameter ^ which roughly shows the fraction 
of segregated states lying in a small vicinity a near ng = 0 
and ng = 1 , 


( 21 ) 

We select here a = 0.1. In Fig. |j the quantity ^ is 
plotted as a function of 7 for various relative strengths of 
residual mixing, panel (a), and for different half-occupied 
systems, panel (b). It follows from the graphs of Fig. ^ 
and § that the segregation starts at 7 ~ D, and occurs 
gradually as 7 dominates the residual mixing V , resulting 
in a peculiar phase transition. The theory of this general 
phenomenon seen earlier in numerical simulations |4^ , 
Q was developed in ||^ where an analogy to the Dicke 
coherent state in optics ||^ was pointed out. A similar 
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FIG. 2: Trajectories of 70 many-body states of the system 
of 4 particles on 8 single-particle levels are shown in a com¬ 
plex plane as a function of the increasing decay width 7 of the 
upper single-particle state. The residual two-body interaction 
matrix elements are selected randomly from the Gaussian dis¬ 
tribution with zero mean and with the variance of one energy 
unit. The single-particle energies are equidistant with a spac¬ 
ing Ae = 0.5 energy units. 


effect, with possible implications for quantum computers, 
is known in the context of quantum measurement theory, 
where coupling and decay created by the measurement 
mechanism can result in a separation of dynamics into 
decoupled subspaces Q. 

Formally the phase transition follows from the factor¬ 
ized structure of the continuum coupling in Eq. (|^). The 
rank r of the matrix W is equal to the number of open 
channels (the dimension of the first group of states in 
the example of Fig. ||) that is smaller than the total 
dimension A/” of the intrinsic space. This matrix has r 
nonzero eigenvalues whereas the remaining Af — r eigen¬ 
values are equal to zero. In a more specific language 
[p^ the independent particle basis |$) of this ex¬ 
ample is the “doorway” basis for the coupling to and 
through the continuum. In the limit of strong contin¬ 
uum coupling, the part W dominates the dynamics and 
aligns the eigenstates along its eigenvectors. Similarly 
to the coupling of two-level atoms through their com¬ 
mon radiation field in the Dicke superradiance, here the 
intrinsic states are coupled through the decay channels. 
At values of 7 exceeding the level spacing (the regime of 
overlapping resonances), this coupling becomes strong, 
and the continuum accomplishes the self-organization of 
intrinsic structure . The real part V of the interaction 
is needed only to establish the correct level density on 
the real axis. The random character of V in the above 
example does not prevent the system from the ordering 
by the continuum. 


FIG. 3: The evolution of the effective occupation number ns, 
Eq. (21), as a function of 7 for the 13 lowest (selected at 
7 = 0 ) many-body states. 


B. Two unstable single-particle levels 

A slightly more complex example is shown in Fig. 

Here we again have eight equidistant non-degenerate 
single-particle orbitals Cj/. Two of them, with different 
quantum numbers, have nonzero widths equal to 7 and 
7 / 2 . With the random interaction turned on, the evo¬ 
lution of the 70 eigenvalues with increasing 7 separates 
four groups of the eigenstates with different decay rates, 
fast, slow and two intermediate. Short-lived states, type 
(a) in Fig. include configurations with both unsta¬ 
ble orbitals being occupied, thus at large 7 their width 
is roughly F^®) ft! 87 / 2 . The number of such conhgura- 
tions, 6!/(2!4!) = 15, is equal to the number of long-lived 
states [type (d)]. The corresponding quasi-stable conhg- 
urations with Fi'^^ ft 0 at 7 ^ 00 do not have a noticeable 
admixture of unstable orbitals so that all four particles 
are distributed over six stable orbitals. Finally, the in¬ 
termediate lifetimes F^^^ ft 7 and ~ 7/2 correspond 
to cases where only one of the unstable orbitals is filled. 
The number of such cases is 6!/(3!)^ = 20 . 


C. Kramers degeneracy 

The manifestations of symmetries and symmetry 
breaking in open systems are particularly interesting and 
important for the shell model formalism. In the pres¬ 
ence of global symmetry, such as rotational invariance, 
the entire space that includes both internal and external 
states, is symmetric. In this case degeneracies of states 
are possible and there is no mixing between the classes 
of states with different exact quantum numbers. A non¬ 
trivial situation occurs when symmetry is violated in ex¬ 
ternal space or in the internal-external coupling, while 
the intrinsic system is still symmetric. Then the degen- 
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FIG. 4: Fraction of segregated states as a function of 7 for the 
same system as in Figs. 3 and 4. For the upper panel (a) the 
single-particle level spacing is varied from the degenerate case, 
Ae = 0, to Ae = 100, a point where the residual interaction 
V can be completely ignored. The lower panel (b) shows the 
“condensed” fraction of many-body states for a degenerate 
single-particle spectrum and various system sizes. In all cases 
only one single-particle level undergoes decay with the width 
7- 


eracy is lifted and in the effective Hamiltonian, in general, 
both non-Hermitian part W and Hermitian part A are 
no longer invariant. Here we consider a schematic case 
when symmetry is violated only in the internal-external 
coupling. 

In systems invariant under time reversal, the states of 
an odd number of fermions are at least double degener¬ 
ate (Kramers degeneracy). Let our eight single-particle 
orbitals form four double-degenerate pairs, as it happens 
in the body-fixed frame of a nucleus with static mul¬ 
tipole deformation where the degeneracy in an axially- 
symmetric case is connected with the sign of the mag¬ 
netic quantum number, ±to. This situation is shown in 
Fig. ||(a) for the case of 3 particles and random selection 
of time-reversal invariant two-body residual interaction. 
All curves in this figure correspond to the evolution of 


-100 -80 -60 -40 -20 0 20 40 so 80 100 



Re(£) 


FIG. 5: Trajectories of 70 many-body states of the system of 
4 particles on 8 single-particle levels are shown in a complex 
plane as a function of the increasing decay width, 7 and 7 / 2 , 
of the two upper single-particle states. The single-particle en¬ 
ergies and the two-body interaction matrix elements are the 
same as in Fig. 2. Arrows (a), (b) and (c) indicate the end¬ 
points at large 7 for the corresponding types of states: short¬ 
lived (a), and two intermediate types (b) and (c), see text. 
The remaining states of the fourth class (d) are quasistable. 


double-degenerate many-body states, since two decaying 
Kramers-degenerate single-particle levels are assumed to 
have the same width, and thus the symmetry is not vio¬ 
lated by decay. 

The case (b) illustrates the situation when time- 
reversal invariance is slightly distorted in the decay chan¬ 
nel, for instance by an external magnetic field, so that 
one of the Kramers-degenerate levels has a 10% larger 
width. The long-lived and short-lived states are the least 
affected ones, in Fig. ^(b) the splitting of their degen¬ 
eracy is hard to resolve. Indeed, for both of these cases 
the decay dynamics either lock a particle pair on the de¬ 
caying time-conjugate orbitals or make them both unoc¬ 
cupied. This restriction of motion allows the system to 
retain the quasi-invariance. The remaining states with 
the intermediate lifetime involve superpositions with one 
particle being on either of the two decaying single-particle 
states. Such a superposition, generally, is no longer time- 
reversal invariant and the Kramers degeneracy is broken. 
Finally, in the limit of strong decay, residual two-body 
interactions become less effective in mixing and then real 
parts of complex energies £ regain the degeneracy. These 
arguments, however, no longer hold true in the extreme 
limit of the violated time-reversal invariance when one 
of the time-conjugate orbitals becomes stable while the 
partner decays. Fig. ||(c). Then the Kramers degeneracy 
is broken and, since there is no intermediate groups, the 
states join long-lived or short-lived families, which leads 
to 35(long)-|-21(short)=56 states. 








































inally unstable states in the presence of the coupling 
through continuum. We consider two single-particle lev¬ 
els, let say Si /2 and pi /2 orbitals, in the continuum so 
that their energies e(s) and e(p) are positive if the con¬ 
tinuum threshold is put at zero energy. In the three-body 
Borromean model for ^^Li with the inert core of ®Li and 
particle-unstable ^°Li, the two active states are those for 
a pair of halo neutrons, ei = 2e(p), e 2 = 2e(s). They 
are quasistationary, and their decay amplitudes Ai ^2 for 
the only open channel, characterized by the core nucleus 
in the ground state and the neutron pair in the state 
= O'*" in the continuum, can be found from a single¬ 
particle picture. At this point the exact form of the en¬ 
ergy dependence is not fixed, except for the fact that 
when the total energy approaches zero, the decay be¬ 
comes forbidden so that, as in Eq. (p^, the amplitudes 
Ai ,2 contain the step function Q{E). A special case with 
one initial non-zero width was presented in [ ^ . The 
problem of two-body decay, especially relevant for Bor¬ 
romean systems, was recently approached with the use 
of different shell model formalisms in Refs. |^, |^. Two- 
proton radioactivity is another example requiring a 
similar consideration. 

According to Sect. II, the effective non-Hermitian 
Hamiltonian in this 2x2 space is 


FIG. 6 : Dynamics of many-body eigenstates of the sys¬ 
tem containing 4 double-degenerate single-particle states and 
N — 3 particles. Panel (a) corresponds to the upper two de¬ 
generate single-particle levels having the same width 7 . The 
Kramers degeneracy is preserved in this case. Panels (b) and 
(c) correspond to the situations where the invariance is bro¬ 
ken by the decay. In (b) the upper pair of initially degenerate 
time-conjugate levels have widths 7 and I.I 7 , respectively. 
The case (c) presents the maximum symmetry violation as 
the time-conjugated levels have widths 7 and 0 . 


V. DYNAMICS OF TWO STATES COUPLED 
TO A COMMON DECAY CHANNEL 


A. Appearing of binding 


f ei - § 7 i u- f ^1^2 \ 

\ U - IA1A2 £2 - |72 / ' 


( 22 ) 


Here Vi 2 = u is the real mixing matrix element, 71^2 = 
Al 2 , and the amplitudes Ai _2 are also real. One should 
be careful with the phases. For a pure internal interac¬ 
tion, the sign of the mixing matrix element V 12 is irrel¬ 
evant, it always can be changed by the redefinition of 
the phase of one of the states, 1 or 2. But with cou¬ 
pling to continuum this change must be accompanied by 
the corresponding phase change in the decay amplitude; 
therefore we cannot simply put Ai _2 = y' 7 i, 2 - 


Here we consider a model that shows how the attrac¬ 
tive real interaction works generating the binding of orig- 


Formal diagonalization of the effective Hamiltonian 
gives the complex energies of the quasistationary states 


J 


"-=2 


ei -l-e2 - -(71 -1-72) 


± 2 - £2) + - -(71 -b 72) - i[(ei - e2)(7i - 72) + ivAiA2] 


1/2 


(23) 


r 


Let us check some particular cases. 


the standard two-level repulsion, 

I 


£± = E± = 


2 L 


ei -b £2 ± \/(£1 - £2)^ + 4z;2 


(i) For the case of stable states, Ai ^2 = 0, we come to 


(24) 
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The lower level reaches zero energy under the condition while from the imaginary part we obtain 


u" = £ 162 . (25) 

(ii) Consider the case with no intrinsic mixing, v = 
0, and two degenerate resonances, ei = £2 = £. Then 
the Hamiltonian consists of the unit matrix e and the 
matrix W of a special factorized type (rank r = 1) so 
that the correct linear combinations are the eigenvectors 
of W ; one of them, the analog of the Dicke coherent state, 
should accumulate the total width, and the second should 
be stable. Indeed, Eq. ( p^ gives in this case 

% % 

£± = £ - -(71 + 72 ) ± -(71 + 72 ) 


r S = e, r = 0, 

\ E = £ - §r, r = 7i + 72 


(26) 


p _ E(7 i + 72) - 71 £2 - 72£i + 2vAiA2 , . 

2 E-£i-£2 


The coupled Eqs. (3^) and (3^) determine E and f. 
For an arbitrary energy dependence of the amplitudes 
Ai 2 {E) that, in order to find the quasistationary states, 
are to be taken in these equations at E = E, this is 
still an implicit solution; even the number of roots can 
change. 

For a sufficiently strong interaction u, the repulsion of 
real energies can bring the lower eigenvalue E_ to zero (a 
threshold value). Then both amplitu des H i ,2 disappear 
together with the eigenwidth r_, Eq. (31b). This means 
that the lowest quasistationary state becomes bound un¬ 
der the same condition (psi). If the mixing increases fur¬ 
ther, the binding energy of the lower state is going down, 


One of such situations with a bound state in the contin¬ 
uum was discussed in a nd explained in terms of the 
effective Hamiltonian in ||15|| . 

(iii) In contrast to the avoided crossing (^) of stable 
levels, the coincidence of two complex eigenvalues is pos¬ 
sible. It requires that two conditions be fulfilled, 


(£1 - £ 2 )^ = 7172 , 

(27) 

(71 - 72 )^ = 16u^. 

(28) 


The coinciding complex energies £±, Eq. (p^), evenly 
divide the trace of the Hamiltonian. In the energy- 
independent two-level Hamiltonian the condition for 
crossing is 


(£1 — £ 2 ) (71 — 72 ) + 4 vAiA2 = 0. (29) 

Which leads to crossing of either real energies £’+ = E- 
or widths r+ = r_ depending if the sign of 

W = (£i - £2)^-k 4 i;^ - i(7i-f 72)^ ( 30 ) 

is negative or positive, respectively. Clearly, in the 
energy-dependent case same remain true for the crossing 
of energies, other conditions, such as for crossing of the 
widths change because generally 7t(E+) ^ 7f(£’_). The 
crossing and anti-crossing of unstable levels were theoret¬ 
ically discussed also in Ref. and experimentally 

studied with microwave cavities PO!] . 

The secular equation for the eigenvalues can be also 
written in a form explicitly separating the real, E, and 
imaginary, T, parts of complex roots (here we again re¬ 
store the tilde sign in order to distinguish the roots from 
the running energy value E). The real part of this equa¬ 
tion gives 

- E(£i -k £ 2 ) - ^(f - 7i - 72 ) + £i £2 -v'^ =0, (31a) 


E_ 


- £i £2 
£1 + £2 


(32) 


As was mentioned in for a similar model with only 
one non-vanishing 7 , this is a prototype of the dynamics 
leading to the binding of nuclei as ^^Li where the residual 
interaction among the valence neutrons is, as we have 
assumed in this Section, of pairing type. 

The higher level, in the point of bifurcation (p^), has 
the energy 



£1 + £2 + \J (£1 + £2)^ + r+(r+ 


71 - 72 )], 
(33) 


where 71,2 are to be taken at energy E = E+. If the 
effective Hamiltonian were energy-independent, both the 
real and imaginary parts of its trace would be separately 
preserved by the complex orthogonal transformation to 
the eigenvectors. This means that we would always have 


r+ -k r_ = tr IE = 71 -k 72 


(34) 


and 

E^ -k E_ = tr £ = £1 -k £ 2 . (35) 

At the bifurcation point, = r_ = 0 , we would have 
r+ = 71+72, A+ = £i+£2, (36) 

while it follows from Eqs. (|^) and ( fHb| ) that 

[A^{E+)^ + A2{E+)^]^ 


r+(E+) = 


£1 + £2 


and 


r-i-(^^-i-) - ii{E+) - i2{E+) 
[A,{E+)^-A2{E+)^f 
£l + £2 


(37) 


(38) 
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in contradiction to the first part of Eq. d^). The 
trace violation occurs because the imaginary parts have 
their own energy behavior with compulsory zero-energy 
thresholds. When the levels are repelled by the mixing 
interaction, their widths are changed by the dynamics 
outside the 2x2 matrix. But the trace is preserved and 
Eq. d^ ) is fulfilled if 


where the higher root £+ is defined by Eqs. (^ ) and d^). 
At low energy, if ^ 0, the behavior of the scattering cross 
section, as well as photonuclear processes, is determined 
by the actual energy dependence of decay amplitudes. 

C. Solutions with energy-dependent widths 


_ jEi 

A2{E) Y ^2 ’ 


(39) 


so that in the entire energy range of interest the two 
partial widths grow proportionally, an interesting excep¬ 
tional case. 


B. Scattering cross section 

In this subsection we consider the scattering cross sec¬ 
tion for the case of two intrinsic states coupled to one 
open channel. Although this cannot be observed with 
the scattering of a neutron pair, the result is relevant 
for the excitation processes of a Borromean system. The 
elastic cross section in the s-wave for a relative momen¬ 
tum k oc '/E is 


a{E) = ^\SiE)-lf, (40) 

where the scattering matrix is defined by Eqs. (i) and d) 
in terms of the effective Hamiltonian Ti.. In our case, Eq. 
(^), neglecting the potential scattering s, the propagator 
can be easily found, and we obtain 


rrfjp\ £'(71-^ 72 ) - 7162 - 72^1 - 2 ?;AiA 2 

=- (£-£+)(£-£_) - ' 

with the poles £± = E± — {i/‘E)T+ given by Eq. (|^), or 
by a pair of coupled equations ( pl[) . One can notice that 
the relative sign of the matrix elements for the direct in¬ 
ternal interaction between the mixed states, v, and for 
their continuum mediated interaction, A 1 A 2 , may con¬ 
siderably change the resulting cros s sec tion. 

In the special case [(ii). Sect. VAj of a pair of de¬ 
generate intrinsic levels with no direct interaction, Eq. 
(p^), the general result (^ reduces formally to the sin¬ 
gle Breit-Wigner resonance on a Dicke coherent state. 


T{E) = 


_7i +72 _ 

E - eE (*/2)(7i + 72 ) 


(42) 


The second root, T = 0, of Eq. (p^ is decoupled from the 
continuum and does not influence the scattering process. 
We have to stress again that the “widths” 717 in general 
depend on running energy E. 

At the bifurcation point (psi), the scattering amplitude 
becomes 


£(71 +72) - (^lv^ + ^ 2 V^)^ 

=- E{E-£^) - 


( 43 ) 


In this and the next subsections we illustrate the dis¬ 
cussed above dynamics of two states coupled to a com¬ 
mon continuum. For all figures here we assume that 
ei = 100 keV and ei = 200 keV for the particle pair 
in p and s states, respectively. For these parameters 
the ground state reaches zero energy, and thus becomes 
bound, at ~ 141 keV by virtue of Eq. (p^. 

The picture with energy-independent widths is not 
consistent with the definition of thresholds. As seen from 
Fig. |(a), the residual interaction pushes the levels apart, 
and the lower state crosses zero energy. However, the 
width r_ of this state, dashed lines in Figs. 0(fa) and 
0(c), is still positive. That would contradict energy con¬ 
servation. For calculations shown by solid lines in Figs, 
^(b) and ^c) we account for the squeezing of the avail¬ 
able phase space volume that forces the decay amplitudes 
to vanish once there is not enough energy for the process 
to take place. Similar to Ref. | |33| , we assume in the low 
energy region the square root energy dependence for the 
s-waves, and ^ E^^^ for the p-wave, 

-/2iE)=aVE, -ii{E) = (3E^/‘^. (44) 


Then the evolution of complex energies as a function of 
the strength v of the residual interaction is consistent 
with the existence of thresholds; at = eie 2 the lower 
state becomes stationary, E- = r_ = 0. The near¬ 
threshold behavior of the width is governed by the s- 
wave component with the infinite slop^ T ~ V£ — EF). 
However, as a becomes smaller. Fig. 0 (c), the singular¬ 
ity is getting confined to a smaller vicinity of threshold, 
to the limit that at an observable scale the behavior is 
dominated by the p-wave. 

Besides the trivial situation, when the width of a par¬ 
ticular state vanishes due to energy conservation, block¬ 
ing of the decay via dynamical mixing at a single point 
corresponding to some strength v is possible. This effect 
of the bound state in the continuum is seen in Fig. 0 (b), 
where a conspiracy of the parameters leads to the vanish¬ 
ing width r_ of the lower st ate a t ene rgy E - still in the 
continuum. Fig. 0 (a). Eqs. (^) and (3^) with r_ = 0 
show that this happens at the interaction strength 


v = A^A2 - —(45) 
7 i - 72 

Here Ai ^2 and 71^2 = A? 9 are to be taken at the energy 
E- found from Eq. (|31cj ). For the model in Fig. 0 (b) 
this happens at a ~ 63 keV. A similar case appears in 
Fig. 0 at c ~ 185 keV, where both states of the system 
become stable only at this particular mixing strength. 
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V (keV) 


FIG. 7: Panels (b) and (c) demonstrate the behavior of F- as 
a function of v for energy-dependent (solid lines) and energy- 
independent (dashed lines) decay amplitudes. Selected pa¬ 
rameters are: Ai = 8.1 (keV)^^^ and A 2 = 12.8 (keV)^^^ 
in the energy-independent case, and a = 15 (keV)^^^ and 
(3 = 0.05 (keV)~^^^ in the energy-dependent case, panel (b); 
Ai = 7.1 (keV)^^^ and A 2 = 3.1 (keV)^^^, dashed line, and 
a = l(keV)^/^ and /3 = 0.05 (keV)“^^^, panel (c). Param¬ 
eters are selected in such a way that at n = 0 the two solid 
and dashed lines agree. The relative phases are such that 
n > 0 and A 1 A 2 > 0 . In panel (a) energies of the two states 
are shown, solid lines, for the case relevant to panel (b) with 
the energy-dependent amplitudes, Eq. (^^, and compared 
to the energies of a non-decaying system, dashed lines. The 
dotted line in all three plots corresponds to the zero value of 
the width or energy. 


The energy dependence of the amplitudes complicates 
the motion of eigenvalues in the complex plane. Interest¬ 
ing features of the level crossing are demonstrated in Figs. 
^ and It should be emphasized that, unlike in a stable 
system or a system with energy-independent parameters, 
here the solutions for £+ and £- involve a diagonalization 
of different matrices. The Hamiltonian matrices differ in 
their imaginary part W. The “interaction” between lev- 


FIG. 8: Possible level crossing in a decaying system. Panels 
(a) and (b) show the energies and widths, respectively, for the 
two-level model of ^'^Li. The upper panel, dashed lines, shows 
also the level repulsion in a closed system. The parameters 
used are the same as in Fig. |(b), except for the opposite 
phase, A 1 A 2 < 0, for n > 0. 


els occurs via common Hermitian part H°+V. Thus for a 
general system it can be expected that bound and weakly 
bound states are still strongly correlated although some 
new features related to small imaginary components ap¬ 
pear, and the usual level repulsion is present only up to 
a spacing of the order of the level width |]^. For states 
deeply in the continuum, however, the correlation must 
rely on the structure of W(E), that represents features 
and symmetries of the continuum. 

Figures ^ and ^ also emphasize the importance of rela¬ 
tive phases for the matrix elements of the internal inter¬ 
action and interaction to the continuum. We can always 
assume u > 0 as was done for our hgures but the change 
in sign of A 1 A 2 leads to a system with very different 
properties, compare Figs. ^(a) and (b) to Fig. ||. 
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FIG. 9: Avoided crossing in a decaying system. See descrip¬ 
tion for Fig. ^ Couplings are modified to q = lO(keV)^^^ and 
/3 = 0.05 (keV)-i/2. 




E (keV) 


FIG. 10: The near-threshold scattering cross section is shown 
for loosely bound systems with v = 180 and 200 keV in 
solid and dashed lines, respectively. Due to similarity be¬ 
tween the curves, the dashed curve is not shown in panel (b). 
Other selected parameters are a = 15 (keV)^^^ and /3 = 0.05 
(keV)~^^^ . Phases are A 1 A 2 > 0 and j4ij42 < 0 for panels 
(a) and (b), respectively. 


D. Cross sections near threshold 


At the critical value of v, Eq. (|^), and in the low 
energy region where the approximation ( p4[ ) can be valid, 
the scattering amplitude (^) is singular, ~ 


T(E) 


aci 

VE{ei + £2 - (i/ 2 )Q:[e 2 /(ei -b £ 2 )] VE} 


(46) 


When the interaction is over-critical, > £i£ 2 , and at 
low energies, E < \E-\ [Eq. (p^], we obtain 


T{E) 


aeiVs 

£+{E-E.y 


(47) 


Therefore the cross section ( |4^ ) has a constant value at 
threshold and behaves at low energies as {E -b |E_|)~^ 
revealing “attraction” to the sub-threshold region [[1^ . 
Similar near-threshold resonance phenomena were dis¬ 
cussed by Persson et al. [p^. The cross sections shown 
in Fig. |^(a) for two over-critical values of the interac¬ 
tion strength reveal a threshold behavior characteristic 
for loosely bound systems that can be mistaken for res¬ 
onances. In this case we had A 1 A 2 > 0 that produces 


only a very broad peak (not shown on Fig. 0 a)) cor¬ 
responding to the upper quasistationary state E+. The 
next Fig. p^(b) shows that in the case of A 1 A 2 < 0 the 
interference of the internal and external interactions re¬ 
sults in a narrow resonance with a very high cross section 
at E = E+ = 340 keV, along with the maximum at zero 
energy (of course, all numerical values characterize only 
the model parameters). 


VI. PAIRING IN THE CONTINUUM 


In this Section we present examples of more realistic 
shell model calculations, where, in general, a large num¬ 
ber of states is involved and further complications arise 
from the conservation of exact quantum numbers in the 
decay as well as from the required self-consistency be¬ 
tween binding energies and thresholds for one and few- 
body decay channels. 

The effective Hamiltonian Ti. implicitly depends on en¬ 
ergy and other quantum numbers that determine if the 
decay is allowed by the conservation laws and what is the 
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near-threshold decay rate. Below we consider examples 
where the intrinsic configuration mixing is generated by 
the pairing interaction only. This leads to the conserva¬ 
tion of all partial seniorities (a number of unpaired par¬ 
ticles Sj on each orbital j). We will also use s = sj 
as total seniority. 


A. Two-level model 

We again start with a simple model of a two-level sys¬ 
tem. We assume here that each level can accommodate 
111 = n 2 = 10 particles and both levels can decay to a fi¬ 
nal state that has fixed Ef = 0 , their decay widths 71 and 
72 are different but have the same, 71 , 2 (^ 1 ) = Q;i, 2 V^, 
energy dependence near threshold. It should be noted 
that introduction of a synchronous decay with 71 = 72 
does not effect the internal dynamics since W in that 
case is proportional to a unit matrix. The corresponding 
single-particle energies are taken as ei = 1 and €2 = 3. 
Intrinsic dynamics in this model are generated by the 
constant pairing V/Iq = G . In Fig. the spectrum of 
states with seniority s = 0 in the system of 8 particles 
is shown as a function of the pairing strength. The at¬ 
tractive pairing interaction pushes down low-lying levels, 
forcing some of them to become bound. The ground state 
becomes bound at G ~ 0.2 . 

Comparison of spectra with and without continuum 
coupling (solid and dashed lines, respectively) shows 
generic features. The bound states are not affected by 
the continuum coupling. The low-lying levels, as com¬ 
pared to highly excited states, are less influenced by the 
presence of continuum. In contrast to usual perturbation 
theory, we see that the ground state and even the first 
excited state once embedded in the continuum become 
attracted to the bulk of other states that increases their 
energy. Such a situation usually leads to an increase of 
the decay Q-value that in turn further increases the decay 
width. 

The following Fig. |l^(a) demonstrates the shift AE of 
the ground state energy as a result of decay for various 
choices of continuum coupling given by parameters oi 
and 02 . Clearly, AE = 0 if there is no configuration 
mixing at G = 0, or once the state becomes bound. The 
complex behavior of the decay width for the ground state 
is shown in Fig. 01 (b); at the critical strength the width 
goes to zero with an infinite slope, ^ y/E. 


B. Realistic pairing model 

As a demonstration of a realistic shell-model calcu¬ 
lation we consider oxygen isotopes in the mass region 
A = 16 to 28. In this study we use a universal sd-shell 
model description with the semi-empirical effective inter¬ 
action (USD) Q. The model space includes three single¬ 
particle orbitals lsi/ 2 , 0 ^ 5/2 and 0 ^ 3/2 with correspond¬ 
ing single-particle energies —3.16354, —3.94780 and 



FIG. 11: The level scheme of s = 0 states in the two-level, 8 - 
particle system as a function of pairing strength. Solid lines 
correspond to the system embedded in the continuum with 
the fixed width values 01 = 0.1 and 02 = 5 . These curves 
are compared with a non-decaying situation qi = 02 = 0 of 
the usual shell-model, dashed lines. The dotted line at E — 0 
indicates the threshold location. 


1.64658 MeV. The residual interaction is defined in the 
most general form with the aid of a set of 63 reduced two- 
body matrix elements in pair channels with angular mo¬ 
mentum L and isospin t, j4T4)Lt\V\{jiTi, j2T2)Lt) , 

that scale with nuclear mass as (A/18)“°’^. 

Although the full shell model treatment is possible 
for such light systems, here we truncate the shell-model 
space to include only seniority s = 0 and s = 1 states. 
This method, “exact pairing -|- monopole”, is known 
to work well for shell model systems involving only one 
type of nucleons (in the case of the oxygen isotope chain 
only neutrons are involved). The two important ingredi¬ 
ents of nuclear forces are treated exactly by this method: 
the monopole interaction that governs the binding energy 
behavior throughout the mass region, and pairing that 
is responsible for the emergence of the pair condensate, 
renormalization of single-particle properties and collec¬ 
tive pair vibrations. In our exploratory study, the trun¬ 
cation of the large space to the most important states 
is a reasonable approach since certainly the inclusion of 
decay makes the computations more numerically intense. 

In the resulting shell-model description the set of the 
original 30 two-body matrix elements in the isospin t = I 
channel is reduced to 12 most important linear combina¬ 
tions. Six of these are the two-body matrix elements for 
pair scattering in the L = 0 channel describing pairing, 
and the other six correspond to the monopole force in 
the particle-hole channel, 

Vjy = E( 2 ^ + mJ,f)Ll\V\U,f)Ll), (48) 

L=ji0 

where j and j' refer to one of the three single-particle 
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FIG. 12: The upper panel shows the shift in ground state en¬ 
ergy between decaying and non-decaying (usual shell-model) 
systems, AS = E{^) — S(4's.m.) , as a function of pairing 
strength under various assumptions for the decay rates. In 
the lower panel the width of the ground state is plotted. 

levels. 

We assume here that the orbital 0 ^ 3/2 belongs to the 
continuum and therefore its energy has an imaginary 
part. In this model we account for two possible decay 
channels for each initial state |$), a one-body channel, 
c = 1, and a two-body channel, c = 2. The one-body 
decay changes the seniority of the 0 ^ 3/2 orbital by one, 
from 1 to 0 in the decay of an odd-A nucleus and from 0 
to 1 for an even-A nucleus. The two-body decay removes 
two paired particles and thus does not change the senior¬ 
ity. The two channels lead to the lowest energy state of 
allowed seniority in the daughter nucleus, i.e. the pos¬ 
sibility of transition to excited pair-vibrational states is 
ignored. This results in 

^ 3 / 2 {^) = £3/2 ~ 2^3/2 ~ 

(49) 

where we assumed that one and two-body decay param¬ 
eters are related as 73 J 2 = 73 ^ 2 /^ — 73 / 2 ) f^e 


particles are emitted in the d-wave with 1 = 2. 

These assumptions can be reviewed by examination of 
where all three states with a valence particle lo¬ 
cated at one of the single-particle orbitals can be clearly 
identihed as the 5/2+ ground state and 1/2+ and 3/2+ 
excited states. Their energies relative to exactly cor¬ 
respond to the single-particle energies in the USD model. 
Furthermore, experimental evidence indicates that the 
3/2+ state decays via neutron emission with the width 
r(^^O) = 96 keV. This information allows us to deter¬ 
mine our parameter 0 : 3/2 = r(^'^ 0 )/(e 3 / 2 )^/^ = 0.028 
(MeV)“^/^ . Other two states are particle-bound, 71/2 = 
75/2 = 0 . 

Using the complex single-particle energies, the effective 
non-Hermitian Hamiltonian for the many-body system 
is constructed in a regular way. We treat the chain of 
isotopes one by one starting from Therefore for each 
A the properties of the possible daughter systems A — 1 
and A — 2 are known. Since the effective Hamiltonian 
depends on energy, and all threshold energies have to be 
determined self-consistently, we solve this extremely non¬ 
linear problem iteratively. We start from the shell-model 
energies i?s.m. corresponding to a non-decaying system 
with the Hamiltonian El . Then the diagonalization of 
EL{Es.ui.) allows us to determine the next approximation 
to the energies. This cycle is repeated until convergence 
that is usually achieved in less than ten iterations. 

The results of the calculations and comparison with 
known experimental data for the chain of oxygen isotopes 
are shown in Table ^ Despite numerous oversimplifica¬ 
tions related to seniority truncation (some widths in the 
Table vanish only due to the fact that only s = 0 and 
s = 1 states were included), limitations on the configura¬ 
tion mixing and restrictions on possible decay channels 
and final states, the overall agreement observed in Table 
I is quite good. The results where experimental data are 
not available can be considered as predictions. In our 
view, however, the main merit of this calculation is in 
demonstrating the power of the method. 


VII. CONCLUSIONS 

The goal of the present paper is to demonstrate a vari¬ 
ety of results that can be obtained with the use of an ef¬ 
fective non-Hermitian Hamiltonian incorporated into the 
standard framework of the nuclear shell model. Although 
those ideas are known for a long time, right now it seems 
to be an appropriate moment to revive them and convert 
into a working tool for the solution of numerous practical 
problems of nuclear, and supposedly more general many- 
body, theory. In all cases when a many-body quantum 
system of strongly interacting particles is loosely bound, 
the interplay of the continuum and intrinsic structure is 
getting crucial, and the phenomena on the borderline be¬ 
tween the bound states and reaction channels become ex¬ 
ceedingly important. Therefore the formalism that would 
allow for a unified description of interrelated structure 
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A 

J 

E (MeV) 

F (keV) 

Eexp. (MeV) 

Fexp. (keV) 

16 

0 

0.00 

0 

0.00 

0 

17 

5/2 

-3.94 

0 

-4.14 

0 

17 

1/2 

0.78 

0 

0.87 

0 

17 

3/2 

5.59 

96 

5.08 

96 

18 

0 

-12.17 

0 

-12.19 

0 

19 

5/2 

-15.75 

0 

-16.14 

0 

19 

1/2 

1.33 

0 

1.47 

0 

19 

3/2 

5.22 

101 

6.12 

110 

20 

0 

-23.41 

0 

-23.75 

0 

21 

5/2 

-26.67 

0 

-27.55 

0 

21 

1/2 

1.38 

0 



21 

3/2 

4.60 

63 



22 

0 

-33.94 

0 

-34.40 

0 

23 

1/2 

-35.78 

0 

-37.15 

0 

23 

5/2 

2.12 

0 



23 

3/2 

2.57 

13 



24 

0 

-40.54 

0 

-40.85 

0 

25 

3/2 

-39.82 

14 



25 

1/2 

2.37 

0 



25 

5/2 

4.98 

0 



26 

0 

-42.04 

0 



27 

3/2 

-40.29 

339 



27 

1/2 

3.42 

59 



27 

5/2 

6.45 

223 



28 

0 

-41.26 

121 




TABLE I: Seniority s = 0 and 1 states in oxygen isotopes. 
Energies and neutron decay widths are shown. Results are 
compared to the known data. Ground state energies relative 
to the core are given in bold. The rest of the energies are 
excitation energies in a given nucleus. 

and reaction aspects is especially needed, and many at¬ 
tempts in this direction made during recent years clearly 
demonstrate this need. 

We illustrated the richness and nontrivial character of 
physics revealed by the complicated interplay of internal 
and external dynamics using a hierarchy of examples, 
from the simplest ones to less obvious to realistic many- 
body problems. Among the most interesting phenom¬ 
ena emerging here we can mention the redistribution of 
the widths, similar to the Dicke superradiance in optics. 


and the segregation of direct processes from those going 
through the compound nucleus stage; interference of the 
“normal” intrinsic residual interaction and the interac¬ 
tion mediated by the excursion into open decay channels; 
dynamics of the poles in the complex energy plane with 
unusual crossings and anticrossings; emergence of bound 
states from the continuum; typical behavior of the reac¬ 
tion cross sections in the presence of loosely bound states. 
The correct account for the threshold singularities of the 
amplitudes of the processes at low energies was an indis¬ 
pensable part of the entire formalism. Finally, we have 
shown how realistic problems of nuclear structure can be 
solved with the aid of this method. In particular, the 
hybrid of the exact solution for the pairing interaction 
with the interaction through the continuum seems to be 
a promising instrument for future development of theory. 

Certainly, the practical implementation of the method 
may be more complicated than in the standard shell 
model with bound states only. The self-consistency prob¬ 
lems of two types, namely (i) a regular solution for the 
complex energies of quasistationary states governed by 
the energy-dependent Hamiltonian and (ii) the consistent 
determination of bound state energies, open channels and 
reaction thresholds for a chain of nuclides connected by 
those channels, may require new computational efforts. 

The main theoretical problem that was not discussed 
above is related to the residual interaction necessary for 
the very formulation of the shell model problem in the 
presence of the continuum. In principle, the effective 
interaction should be energy-dependent and complex; it 
has to be consistent with the rest of the shell model input, 
including the amplitudes of the coupling to closed and 
open channels. This is a serious challenge for the future 
that requires a new insight into the whole physics on the 
borderline between structure and reactions. 
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